COMPLEXITY OF RANDOM SMOOTH FUNCTIONS OF MANY 

VARIABLES 



ANTONIO AUFFINGER AND GERARD BEN AROUS 

Abstract. We analyze the landscape of general smooth Gaussian functions on the 
sphere in dimension N, when N is large. We give an explicit formula for the asymp- 
totic complexity of the mean number of critical points of finite and diverging index at 
any level of energy and for the mean Euler characteristic of level sets. We then find two 
possible scenarios for the bottom energy landscape, one that has a layered structure of 
critical values and a strong correlation between indexes and critical values and another 
where even at energy levels below the limiting ground state energy the mean number of 
local minima is exponentially large. These two scenarios should correspond to the dis- 
tinction between one-step replica symmetry breaking and full replica-symmetric breaking 
of the physics literature on spin glasses. In the former, we find a new way to derive the 
asymptotic complexity function as a function of the IRSB Parisi functional. 



1. Introduction 

This work deals with the complexity of smooth functions of many variables. The ques- 
tions addressed in this paper can be phrased as: What does a random Morse function 
look like on a high-dimensional manifold? How many critical values of given index, or 
below a given level? What can be said about the topology of its level sets? We study here 
general smooth Gaussian functions on the sphere in dimension N, when N is large. We 
investigate the number of critical points of given index in level sets below a given value, as 
well as the topology of the level sets through their mean Euler characteristic. Our main 
result is that these functions have an exponentially large number of critical points of given 
index, and that the Euler characteristic of the level sets have a very interesting oscillatory 
behavior. Moreover we find an invariant to distinguish between two very different classes 
of complexity for these functions. These two classes should correspond to the distinction 
between one-step replica symmetry breaking and full replica-symmetric breaking of the 
physics literature on spin glasses. Indeed the general random Gaussian smooth functions 
on the sphere correspond exactly to the Hamiltonians of an important class of models of 
statistical physics of disordered media, i.e. mixed spherical spin-glasses. These mean-field 
models, as well as other spin glass models, are well-known to be very challenging to ana- 
lyze. It is believed (see |CL04) and the references therein) that a subset of the spherical 
models that we study here share the same interesting static and dynamical behavior as 
the famous Sherrington-Kirkpatrick model at low temperature. As part of our study, we 
give further evidence for this claim and conjecture its domain of validity. 

We start by calculating the averaged complexity of these functions i.e. the exponential 
rate function of the mean number of critical points of finite and diverging index at any 
level of energy. This initial computation uses the method developed in [ABGlOj . where 
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this study was initiated for particular covariance functions that appear as Hamiltonians 
of the pure spherical p-spin models. 

In the general case, our first result is that the complexity of critical points of finite 
index can be decomposed into two pieces: only one of them present in the pure case (see 
Figure 1). This difference allows us to separate the models of Gaussian smooth functions 
on the sphere in two classes : one where the bottom landscape is qualitatively similar to 
the pure p-spin models and another which should correspond to the case of Full replica 
symmetry breaking where in particular even at energy levels below the limiting ground 
state energy, the mean number of local minima is exponentially large. 

In the former case, that we call pure-like region, we prove a strong correlation between 
critical values and their indexes. There exist energy thresholds —E^ such that below —Ek 
with probability going to one it is only possible to find critical points of index less than 
k. In the latter case, called full mixture region this layered structure is not present and 
there is no difference between the complexity of critical points of finite index k for any k 
(not diverging with A^). However, in both cases the complexity of critical points of any 
index is different from the pure case as we increase the level of energy. In particular, 
coexistence of local minima and local maxima is perfectly possible and the mean number 
of critical points of any finite index agree in a full neighborhood of energies around their 
"most typical" energy (see Theorem 2.1). 

The understanding of the landscape of these Hamiltonians might prove useful for the 
study of both static and dynamical questions of these models. First, the layered structure 
described above may shed a light on the metastability of Langevin dynamics (in longer 
time scales than those studied in [BG97]). Second, it may provide an insight on the most 
important statics open question: to understand at any temperature the structure of the 
(random) Gibbs measure associated to these Hamiltonians. In this direction, a major 
breakthrough was done by Talagrand |Tal06a| based on the remarkable work of Guerra 
|Gue03| by computing the free energy at any temperature under a convexity assumption. 
The free energy is given as the infimum of the Parisi functional over the space of probability 
measures on [0, 1] (see ( |3.2[ )). Understanding the minimizer of this functional (uniqueness, 
for instance, is only known under certain conditions), called a Parisi measure, is also a 
major challenge. In the spherical pure p-spin it was shown in Proposition 2.2 of |Tal06a) 
that the model has a one step replica symmetry breaking (1-RSB) at low temperature, i.e 
the Parisi measure at low temperature is atomic with two atoms. For a mixed spin model, 
as far as we know, the structure of the Parisi measure remains an open question. 

We also show that in the pure-like region, the 1-RSB picture is consistent with the 
complexity picture without any convexity assumption. Precisely, the complexity function 
of local minima can be characterized near its zero as a function of the Parisi Functional 
minimized over two-atomic measures and vice- versa. This is the content of Theorem 13.81 
Furthermore, we show that concentration of the number of local minima implies 1-RSB at 
zero temperature for any mixture in the pure-like region. In the full mixture region, the 1- 
RSB Parisi functional does not describe the complexity function near its zero. In Theorem 



3.1 we show that the ground state energy converges almost surely to a constant. In the 



full mixture region, this constant has positive complexity. This immediately implies that 
it is not possible to have concentration for the number of critical points around its mean 
and one may expect that the averaged complexity is strictly bigger than the quenched 
complexity (see the discussion after Theorem 3.8 ) 

Our picture is consistent and generalizes the one proposed by physicists. In |CL04] . 
it is claimed that a 2 -\- p spherical spin glass model with p > 4, at low temperature is 
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either 1-RSB or its Parisi measure has an absolute continuous part (a Full RSB or a 1- 
Full RSB) depending on how much weight is assigned to the 2-spin model. The regions 
pure-like and full mixture seem to numerically agree and to extend (since we do not need 
the 2 spin component) the one proposed by |CL04| . We find this a remarkable fact : some 
mixtures of the spherical model are expected to have the same Gibbs structure as the 
Sherrington-Kirkpatrick Model on the hypercube. We conjecture that Full RSB holds in 
the full mixture region. Intuitively, since we prove that the average number of critical 
points at the ground state energy is exponentially large with N, the Gibbs measure at 
low temperature has plenty of candidates to sample from. However, our techniques are 
still worlds away to prove this fact. In particular we still do not know what is the typical 
overlap of two critical points. 

Back to the topology of level sets, we show that the total number of critical points 
at a given level of energy is asymptotically equal to the number of critical points of a 
particular index (that depends on the level of energy). Loosely speaking, at lower levels, 
local minima dominate. In a certain threshold energy window {—NEoo, NEoo) critical 
points of diverging index give the main contribution to the total complexity. Above NEoo , 
the total complexity is equal to the complexity of local maxima. 

This phenomena is related to the asymptotic mean Euler (or Euler-Poincare) character- 
istic of the set of points below a certain level of energy. We show that in absolute value the 
mean Euler characteristic of these level sets is asymptotically equal to the number of total 
critical points at that level. It is therefore exponentially large for most energies. Moreover, 
we prove that is positive outside the window {—NEoo, NEoq) but oscilates 0{N) times 
from positive to negative (exponentially large) values inside {—NEoo, NE^o)- We find this 
picture very interesting but quite hard to visualize. 

The paper is organized as follows. In section [2j we define the model and we state the 
main results about the complexity of critical points. In section [3j we state a few relations 
between the structure of the Parisi measure, the global minima of the Hamiltonian and 
the complexity of critical points. We also define the regions pure-like and full-mixture. 
Next, in section |4j we state our results about the Euler characteristic of level sets. In 
section [5] we prove all Theorems about the complexity function. Their proofs follow the 
same strategy of |ABC10j . Namely, they will follow from an exact formula for the mean 
number of critical points of index k that translates the problem to a Random Matrix 
Theory question. This formula is more envolving than the pure case since in a mixture 
the Hessian matrix gains an independent Gaussian component on the diagonal. This leads 
to different variational principle that we analyze. To end, in sections 6 and [7] we prove the 
results of section [3] while in section |8] we prove the results of section 4 
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2. Complexity and Energy Landscape 

The state space of the spherical spin- glass model is S^~^{^/N) C M^, the Euclidean 
sphere of radius ^/N. A configuration <t is a vector of satisfying the constraint 

i=l 

The Hamiltonian of the pure p- spin model is the random function defined on S"^^^ ( ViV) 

by 

1 ^ 

HnA^) = jv(P-i)/2 E Jh,...,i,<^n---^i„ (T = {au...,aN)eS''-\^), {2.2) 

ii,...,ip=l 

where Ji^^...^ij, are independent centered standard Gaussian random variables. 

Equivalently, H^^p is the centered Gaussian process on the sphere S^~^{^/N) whose 
covariance is given by 

N 

E[Hn,p{^)Hn,p{<^')] = Afi-P(^a,a^)'' = NR{a,a')P, (2.3) 

1=1 

where the normalised inner product R{(t, cr') = ^{cr, a') = ^ YliLi '^i'^'i usually called 
the overlap of the configurations cr and cr' . 

The study of the landscape of such Hamiltonians was done in |ABC10] via a Random 
Matrix Theory. In this paper, we will consider the analogous analysis for a mixed p-spin 
model, i.e. linear combinations, of these Hamiltonians. At a first sight, this question 
appears just as a simple generalization, however, as mentioned in the introduction we find 
a richer structure when considering mixed Hamiltonians instead of pure p-spins. 

We now define a mixture of p-spins (or a mixed spin). Given a sequence (3 = (/3p)pgN,p>2 
of positive real numbers such that 

oo 

J2'^P(3p<oo, (2.4) 

p=2 

let 

oo 

^JV,/3(cr) = ^/3pi/^,p(o-), (2.5) 

p=2 

where for any pair of values p ^ p' the Hamiltonians Hn,p, H^y are independent. Con- 
dition (2.4) is more than enough to guarantee that the above sum is a.s. finite and the 
Hamiltonian Hj^^i^ is a.s. smooth (see Theorem 11.3.1 of [AT07| ) . 
In this case, we have that 

oo 

1K[HnM<^)HnA^')] =NY,fp(R{a,a')Y = Nu{R{a,a')), (2.6) 

where 

oo 

Kt) = E/3'*'- (2-7) 

p=2 

A word of comment is needed here. By Schoenberg's theorem [Sch42| . if L'{R{a,cr')) is 
a positive-definite function for all N and all <t, cr' G S^~^{\nV) then u can be written as 
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a linear sum as in (2.7). This shows that we are exhausting all possible covariances given 
as ( [2^ . 

From now on, we call the function v a mixture and we note that v is smooth with 

oo 

v'{\) = 0, ^u">0 = Y,Pl = 1. (2.8) 

p=2 



If we consider the random variable X that assign probability j3p to the integer p, then its 
probability measure is given by = ^ f^p^p and 

W.X = v' and = \3.i X = u" + u' - u'"^ . (2.9) 

A mixture is pure if and only if a = 0. Furthermore, note that v" > v' with equality only 
in the pure case with p = 2. The parameters v' ^ v" and o? will be fundamental in our 
analysis. 

We now introduce the complexity of spherical spin glasses as in [ABClOj . For any Borel 
set i? C M and any integer < A: < A/^, we consider the (random) number Crt7V fc(i?) of 
critical values of the Hamiltonian -f^Ar,/? in the set NB = {Nx : x £ B} with index equal 
to k, 

CTtN,k{B)= Yl l{HN,(3{a)€NB}l{i{V^HNM^)) = k}. (2.10) 

<t:VHjv,0((t)=O 



Here V, are the gradient and the Hessian restricted to S^~^{VN), and iiV"^ H n ^^{ct)) 
is the number of negative eigenvalues of the Hessian V^i^Tv.pi called the index of the Hessian 
at cr. We will also consider the total number Cii^{B) of critical values of the Hamiltonian 
i^TV,/? in the set NB (whatever their index) 

CiiN{B)= Yl 1{Hn,/3{<t) e NB}. (2.11) 

T:V/f]v,/3(<T)=0 

Our first results will give exact and asymptotic formulas for the mean values E CrtAr ^(i?) 
and ECrtAr(i?), when A — )• oo and k, B and u are fixed. We will compute the limits of 
^ logE Crtjv,fc(-B) and logE CrtAr(i3) as A tends to infinity. 

2.1. Complexity functions for critical values of finite index. 

Our first result is the existence and characterization of the asymptotic complexity of 
the mean number of critical points of index A; in a certain level of energy. 

Theorem 2.1. For any fixed integer k > 0, there exist a continuous function O^^uiu), 
called the k-complexity function, explicitly given in (5.10), such that, for any open set 
SCR, 

lim 4logECrt7v,fc(S) = sup0fc,^(n). (2.12) 



We decide to postpone to section 5^ the explicit expression of the k- complexity func- 
tions 9k,u{u)- However, we describe some important properties of these functions (see 
Figure 3) in the proposition below. We first fix three important thresholds that depend 
on u. Let 

v' + 

and 



E'o. = -r±-i,. E^= . , ^ (2.13) 
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Figure 1. /^-complexity functions Ok^y{u) for —6 < u < —1, k = 1,2,3,5 
in the case where u is pure-hke, i.e. 9k,i^{—Eoo) > 0. The dashed Une 
is the continuation of the parabola that describes Ok,v{u) in the interval 
[— E'oojOo) where they all agree. 



E» = ^ . (2.14) 

Note that 

E^<E'^< E^. (2.15) 
Furthermore, E'^ = E^o if and only if Eq^ = E'^ if and only if u" + v' — i/'^ = 0, that is, 
any equality in ( 2.15| ) implies a triple equality. It occurs if and only if the mixture is a 
pure p-spin (see (2.9)). 

Proposition 2.2. For any mixture v and any k > 0, the k-complexity functions 9k,u{u) 
satisfy the following: 

(1) 9k^^{u) is continuous on M and difjerentiable on M \ {—Eoo}. 

(2) 0k,uiu) is strictly increasing on (— oo, —E'^) and strictly decreasing on {—E'^,oo). 
Its unique maxima is independent of k and equal to 

:= 9kA-E'^) = ^ l°g ^ - Jj!^ > 0- (2-16) 

(3) 0k,v{u) has exactly two distinct zeros. The largest zero is given by —E^ and 
therefore is independent ofk. 

(4) For any k,k' £ N with k < k' , 6k,u{u) < Ok\y{u) for all u G (— oo, —Eoo). 

(5) For any k,k' gN with k < k' , 9k,v{u) = 9k'^v{u) for all u E [—Eoo-,oo). 



Immediately from Theorem 2.1 and Proposition |2.2| we obtain: 
Corollary 2.3. The mean total number of critical points of index k satisfies 

lim 4logECrt^,fc(M) = S^. (2.17) 
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Furthermore, if B = {—oo,u) with u < —E'^ then 



lim ^logECitN,k{-oo,u) = ek,u{u). (2.18) 



Remark 1. By symmetry, Theorem 2.1 also holds as stated for the random variables 
CvtN,N-i{B), with / > 1 fixed if one replaces 6k y{u) by Ok^y{—u). 



We now use Theorem 2.1 and Proposition |2.2| to describe the bottom landscape of the 
mixed spin glass models. For any integer A; > 0, we introduce = Ek{v) > as the 
unique solution in (E'ooiCo) to (see Figure 3 again) 

ekA-Ek{^)) = ^- (2-19) 

That is, —Ek{i') is the smallest zero of the fc-complexity function. It is important to note 
that, by items (4) and (5) of Proposition 2.2 the sequence {Ek{i'))kGN is non-increasing. 
Its structure is of extreme importance and will be further explored in the next section. At 
this point we have the following consequence of Theorem |2.1[ 



Theorem 2.4. For k > and e > 0, let Ajv,fc(e) to be the event "there is a critical value 
of the Hamiltonian Hjyj^ below the level —N{Ek{u) + e) and with index larger or equal to 
k", that is 

oo 

Mkie) = {X]Crt7v,i(-^fc(z^) -£) > 0} 

i=k 

and BN^i:{e) be the event "there is a critical value of index k of the Hamiltonian Hf^^^ 
above the level —N{E^ — e) ", that is 

BnA^) = {Crt7v,fc((-^^ +e,oo)) > 0} 
Then for all k >0 and e > 0, 

limsup^logP(yljv,fc(e)) < and limsup ^ logP(BAr,fc(e)) < 0. (2.20) 



Theorem 2.4 says that with overwhelming probability all critical values of the Hamilton- 
ian Hj\[ fj of index k are inside the interval [—NEf^, —NE^]. A similar result was derived 
for the pure spin glass models in [ABClOj . However, in the pure case it was shown (Theo- 
rem 2.2 of [ABClOj ) that the probability of finding a critical point of finite index above the 
level —NEoo is asymptotically of order exp(— A^^C). Hence, in the mixture case not only 
the window of possible values for the Hamiltonian for a critical point of finite index has 
changed, but also the probability of being outside of that window is of order exp(— A^C). 

2.2. Complexity function for critical values of diverging index and the total 
number of critical points. 

We end this section by studying the number of critical points with diverging index and 
the total number of critical points (regardless of index). Let k = k{N) be a sequence of 
integers such that as A^ goes to infinity 

k{N) 



N 

Let Sy £ (— \/2, \/2) be defined as solution of 



7e(0,l). (2.21) 



1 /■"'"^ / 

- / \/2-x2dx = 7. (2.22) 



TT 
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The first result in this subsection is the analogue of Theorem 2.1 for critical points of 
diverging index. 



Theorem 2.5. For any sequence k{N) satisfying (2.21), as N goes to infinity 



lim — logECrtiv,fc(Ar)(5) 



1 ly" 1 

sup <^2^og— + 



sup 6^,uiu)- 

y&B 



2u" 



a 



2 



(2i/" 



Remark 2. From Theorem 2.5 one can easily get analogues of Theorems 2.4 and Corollary 



2.3 for the case of critical points with diverging index. Its statements are adapted rewrites 
of the respective results. We leave it to the reader. 

We also provide the complexity for the expected total number of critical values at a level 
of energy. Our next result can be described as follows: the mean number of critical points 
at level u is asymptotically given by the mean number of local minima, local maxima or 
critical points of index k{N) ~ 'y{u)N if n < —E'^,u > E'^, —E'^ < u < E'^, respectively. 
Here, j{u) € (0, 1) is such that = V^^, see (2.22). 

Precisely, define 



{9oAu) if u<-E'^ 
OoA-u) if u>E'^ 



-u'+u" ' 



) = sup^g(o,i) = 6'7(«),i/(n), otherwise. 



Theorem 2.6. The total number of critical points satisfies 

lim ^ logE CrtAr(i?) = sup 9y{u) := Qy{u). 



(2.23) 



(2.24) 



3. The Ground State Energy, Pure-like and Full Mixtures and 1-step 

Replica Symmetry Breaking. 

The goal of this section is to establish relations between the structure of the Parisi mea- 
sure, the global minima of H]\f j3, and the results obtained for the asymptotic complexity 
in last section for different classes of mixed spin glass models. We start by recalling known 
results about the free energy at positive temperature, more precisely the Parisi formula 
as proved by Talagrand [Tal06a]. 

3.1. The Parisi Functional. The partition function of the p-spin spin glass is given by 

ZnAP) = I e-^^^'.''('^)A,v(d^), (3.1) 

where Ajy is the normalized surface probability measure on the sphere S^~^{^/N). Let 
M[0, 1] the space of probability measures on [0,1]. By Theorem 1.1 of [Tal06a] . if v is 
convex, the following limit holds almost surely, 

Foo(/3):= hm i-logZ^,,(/3)= inf F,(/3,p) (3.2) 

N^oo JS peAf [0,1] 
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A formula for Fu{f3,p) is given in (1.11) of |Tal06a| and we reproduce it now. Given a 
probability measure p on [0, 1] consider its distribution function Xp : [0, 1] — >• [0, 1]. Write 
for q £ [0, 1]: 

x{q) = / Xp{s)d s. (3.3) 



q 



Assuming that x{q) = 1 for some q < 1, then 



FuiP, p) = 2 [/^ Xpiqy{q)d q + ^+ log(l - q) ] . (3.4) 



If g = 1, we set Fy{j3,p) = oo. A measure that minimizes the right side of (3.2) is called 
a Parisi Measure. It is believed that the nature (atomic, absolutely continuous) of the 
Parisi Measure is of extreme importance to understand the statics of the spin glass model 
|Tal06aL[Ta06b| . 



It is a difficult task to handle the infinite dimensional variational principle in (3.2). 
However, in some cases, ( |3.2[ ) can be simplified. Let Mi [0,1] be the space of atomic 
probability measures on [0, 1] that have at most 2 atoms. Define 

Fi(/3):= inf F,(/3,p). 

peAfi[o,i] 

Clearly, -Foo(/3) < Fi{[3). When equality holds, that is when i^i(/3) = Foa{f3), we say that 
the model has a 1 step replica symmetry breaking (IRSB) at inverse temperature /3. 
In |Tal06aj , it was shown that the pure p-spin glass model is 1-RSB for /3 sufficient large if 
p is even. It is an open question to determine mixtures and values of /? that IRSB holds. 
Let 

GS^ = ^ inf HnA<^) (3.5) 

be the normalized absolute minima of the Hamiltonian, i.e., the energy of its Ground 
State. A straight forward exercise shows that, for any e > 0, with probability approaching 
one as N goes to infinity, 

- Eo{iy) - e < GS^ < limM -^F^ifS) + e. (3.6) 

/3-s>oo p 

The main question we investigate in this section is whether the lower and upper bounds 



given on (3.6) are optimal, that is, whether is possible to identify the limit ground state 
energy using the partition function and the asymptotic complexity. The question to find 
the Ground State Energy is one of the foundational and most relevant questions in the 
study of spin glass system among probabilists |Tal03| Chapter 1] . The left and right sides 



of (3.6) are quantities that come from different computations, complexity and free energy, 
respectively. So, a priori, there is no reason to expect that these bounds match. We start 
our analysis by the following fact: 



The upper bound in (3.6) is optimal in the sense that 



Theorem 3.1. For any convex covariance function v: 
(1) The following limit exists 



lim iFoo(/3) = /oo G [0,oo). (3.7) 



(2) The ground-state energy GS^ converges almost surely to — 
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We turn to ask the same question about the lower bound. Our approach is simply to 
try to prove that foo = Eq{i') directl}|^ The problem is that handling /oo is not an easy 
task since it comes from a rather complicated object, the infimum of the Parisi functional 
over the space of probability measures on [0, 1]. Instead our approach is to compare Eo^iy) 
to the analogous constant as if the model was 1 RSB at low temperature. We will show 



in Lemma 7.1 that for any mixture the following limit exists 



/i := lim = inf {-(b + u'a+lilog^^)]], (3.8) 

where e is a positive constant depending on u. When /i = f^o we say that the model 
is 1-RSB at zero temperature. Surprisingly, the comparision between /i and Eq will 
heavily depend on the structure of the mixture and on the bottom landscape of Hf^j^. 

3.2. Pure-like mixtures and full mixtures. 

In this subsection we relate the Parisi Functional to the asymptotic complexity of spin 
glasses and derive more precise information about the landscape of H^^f^. We first identify 
the regions of mixtures mentioned in the introduction. We refer the reader to Figure 3. 
Let 

G{u',u") := log^ - ^- ""'^"^ ^ ' = eoA-Eoo). (3.9) 

Definition 3.2. A mixture v is called a pure-like mixture if and only if G(i^', v") > 0. If 
G{i'' , v") < 0, is called a full mixture. When G{v\ v") = 0, is called critical. 

Example 1. One can easily verify that all pure p— spins, v{x) = x^, p > 3 are pure-like 
while the spherical SK model, p = 2, is critical. A picture of these regions is given in 
Figure 3. 

Example 2. Consider the case 

i^{t) = ^it^ + (1 - ^ly (3.10) 

where /i G [0, 1]. Then, if p > 3 then it is possible to show that there exists a < ;Uc(p) < 1 
such that V is pure-like if and only if /i < f^dp)- l^c{p) is given as the unique zero in (0, 1) 
of 

(p2-2p)(l-^)(2(p2_p)_3(p2_2p)^+(p_2)V) 1, r, 2p 



+ log 



l+P 



p-\-2fi- p^ 



2((p2_p)(i_^) + 2^)(p + 2/i-p/i)2 2 

see Figure 2. Remarkably, p = 3 in ( |3.10 ) is the only case where the mixture is a pure-like 
mixture for all values of t. 

Our first statement concerning pure-like mixtures is the following result about the bot- 
tom landscape. Let 



2vWu" + -^4z^'V'2 - {u" + v'){2[u" -ly' + v>'^) - {u" + u' - u'"^) log ^ 



It follows directly from the definition of pure -like and (2.19) that: 

Proposition 3.3. If v is a pure-like mixture then the sequence Ej.[u) is strictly decreasing 
and Ek{u) converges to E^ as k goes to infinity. 



^One could argue that if we prove a concentration result for the number of local minima then Eq{v) = 
/oo = limGSjv a.s.. Unfortunately, we still do not know how to control the second moment of Crto,^, since 



we could not derive a manageable formula like (5.2 1 nor we could remove the expectation in Theorem 2.1 
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Figure 2. Function G{v' , v") in the case v = /it^ + (1 — 



Last Theorem combined with Theorem 2.4 says if the mixture v is pure-like then the 



landscape of v at low levels of energy is similar to the pure case. In particular, the 
same interesting layered structure for the lowest critical values of the Hamiltonian Hf^j^ 
holds. Namely, the lowest critical values above the ground state energy are (with an 
overwhelming probability) only local minima, this being true up to the value —NEi{u), 
and that in a layer above, {—NEi{v), —NE2{iy)), one finds only critical values with index 
(local minima) or saddle point with index 1, and above this layer one finds only critical 
values with index 0, 1 or 2, etc. 

Using the fact that /i comes from an easier variational principle, Lemma 5.3 of |ABC10) 



shows that miraculously in the pure case with i/(rc) = xP,p even, (3.6) is optimal as indeed 
we have Eq^i/) = f^o = fi- This result extends as: 

Theorem 3.4. If is pure-like or critical then fi = Eq{v). 



Combining Theorems 3.1 and |3.8| we have the following: 



Corollary 3.5. In the case of v pure-like or critical, concentration of Citi\ifi{—oo,u) 
around its mean implies 1-RSB at zero temperature. 

A word of comment is needed here. It is reasonable (although we do not have a proof 
at the moment) to believe that log CrtAr^o(~oo, n) concetrates around its mean. How- 



ever, in Theorem 2.1 we study an "averaged" complexity instead of the possible smaller 



"quenched" complexity: 



lim ^Elog CrtAro(— oo, n). (3.11) 



It is not clear if (and when) both quantities agree. We conjecture that in the pure-like 
region "quenched" is equal to "averaged" and indeed we have 1-RSB at zero temperature. 
We also believe that when i/ is a full mixture the averaged complexity is indeed larger than 
the quenched complexity. These conjectures are supported by physicists |CL04j and by 
the following result that tells us that the complexity of minima can be constructed from 
the 1 RSB Parisi functional at zero temperature and vice-versa. 
For b £ (0, oo) define 

f^{b)= inf |l(6 + ^'a+^(log^))|. (3.12) 
ae(o,oo) {2 \ a / ) 
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and set 



^^'-2 \-xfi{x), x>Cu , . 

9i{x) = { ^ , ^ ^ (3.13) 



^/v'[v' - 1) \-Cufl{Cy), X<Cy. 

Theorem 3.6. If v is 'pure like then for all u < — i?. 



eo,y{u)= min [ub-bfi{b)]. (3.14) 

Moreover, gi{x) is a convex fucntion, strictly convex in (ci,,oo) and if we set for u > Eoo, 
tpiu) = —Oo^ui—u) then ■0 is the Legendre-Fenchel conjugate of gi{x): 

Tpiu) = max ( ux — gi{x) I . (3.15) 

Remark 3. In fact, such duahty is the first sign at zero temperature of an apparently deeper 
connection predicted by physicists |CGG99| ICGPM05] between the Parisi Functional at 
finite temperature and the TAP complexity (see |ABC10j . section 6 for a definition of TAP 
equations). We plan to explore this connection in the future. 

The above connection does not hold if i/ is a full mixture. We end this section by 
analyzing this case. 

If V is either critical or a full mixture it follows from Theorem |2.1| that for any k finite 
the mean number of critical points of index k are asymptotically equal at any possible 
level of energy. In particular, 

Corollary 3.7. // u is either critical or a full mixture then for any k, k' £ N, 

Ek{i^) = Ek'{i^) = Eo{i^). (3.16) 

Furthermore, for any fc E N the probability of finding a critical value of index k below the 
level —N^Eq^v) + e) is exponentially small in N. 

Theorem 3.8. If i' is full mixture, not necessarily convex, then f\ < Eq{v). 



If is a full mixture, last theorem combined with Theorem |3.1| immediately implies 
that —Eq{i') < —fi < GSn for large enough with probability one. Hence, we can not 
remove the expectation from Theorem 2.1 Moreover, 

Corollary 3.9. If i' is a full mixture, then for any u € {—Eo{i'),—foo) the probability 
of having a critical value below u goes to zero while the mean number of local minima is 
exponentially large in N . Namely for such u there exist constants < Ci < C2 such that 
for N sufficiently large 

ECrtjv,o(-oo,n) > e^^\ and p( Crtiv(-oo, u) > e^^M < e~^^\ (3.17) 



4. EULER CHARACTERISTIC OF LEVEL SeTS 

In this section, we investigate the landscape of the Hamiltonian H^^f^ by analyzing the 
mean Euler characteristic of level sets as N goes to infinity. In order to state our results 
we need further notation. The Hermite functions (pj, j £ N, are defined by 

cPj{x) = (2^j!0f)-i/2/i,(x)e-^ , (4.1) 
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Figure 3. Graph of u' x . In blue, the level set G{v' = i.e. the 
case where v is critical. Dotted lines are the possible values of {v'^iy") for 
the mixtures 2 + 6,2 + 10 and 4 + 30. The gray region is outside the domain 
of possible values for {v\ u"). 



where /ij, i G N are Hermite polynomials, 

h,{x) = e^\-j-ye-^\ (4.2) 

In particular, /io(x) = = 2x,h2{x) = 4x^ — 2x. The Hermite functions are or- 

thonormal functions in M with respect to Lebesgue measure. 
We denote by xi^u) the Euler characteristic of a level set 

An := {cr G S^-\VN) : HnA^^) < ^u}. 

x(-) is a topological invariant, integer valued function that is defined for any CW-complex 
as the alternating sum of Betti's numbers |War71| . It is a functional that is invariant 
under homotopies and satisfies 

x{AuB) = x{A)+x{B)-x{AnB), xW = 1 and ^(^tv) = 1 + (-l)"^^' (4.3) 

where B denotes a A^-dimensional unit ball, Sn the A^-dimensional unit sphere and A, 
B are CW-complexes. x(") roughly measures the number of connected components and 
its number of attached cylindrical holes and handles. Since we are only interested on 
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Euler characteristics of level sets of almost surely Morse functions, we use the equivalent 
definition that follows from Morse's theorem (see (AT07[ Theorem 9.3.2]) : 



N-l 



k=0 



The strategy of using Rice's formula to compute Euler characteristics of level sets was 
developed in [ ATOTl ITTA051 ITAOSj and also explored in [AW09J . In fact, in a similar 
fashion, we prove the following proposition: 



Proposition 4.1 

E 



dxdy. 
(4.4) 



Our main result in this section is the asymptotic formula for Ex(^u) and its relation to 
the asymptotic complexity of the total number of critical points (see (2.24)). 



Theorem 4.2. The mean Euler-Poincare characteristic Ex(^„) satisfies the following: 
(1) Ifu<-E'^, 



lim ^logEx(^„) = e^(u). 
(2) // -E'^ <u<0, with u = -E'^ cos uj G (0, vr) 

247r2iV4 /(a;)(sinu;)V^ 



(4.5) 



Nt{uj) + p{uj) 



l + 0(7V-i)). (4.6) 



where 



1 



sin 2bJ — 2ijj) 



1 Svr 

+ ^ + a(w), 



c{N,u) is given in ( 8.22[ ) and f{uj), a{u}) are given in ( |8.23 ). 
(3) Ifu>Owe have Ex(A„) = Ex{A^u) for N even and Ex(A„) = 2 - Ex(A-u) for 
N odd. 



Let us describe in words the landscape picture emerging from Theorem 4.2 Roughly 



speaking. Theorem 4.2 says that the mean Euler Characteristic of A^ is in absolute value 
asymptotically equal to the total number of critical points at level Nu if u < Eq. This 
picture is fairly intuitive and easy to explain in the bottom of the landscape. As we increase 
the energy level u from negative infinity to —E'^, the level set A^ is "essentially" a union 
of disjoint simply connected neighborhoods of local minima. Since these are exponentially 
large and dominate the total number of critical points, the mean Euler characteristic is 
positive and of same size. As we cross the level —E'^, local minima cease to dominate. 
The total number of critical points and the Euler characteristic (in absolute value) is 
given by the critical values of dominant divergent index. The landscape is then hard to 
visualize. By increasing a tiny amount of energy it oscillates from a large positive to a 
large negative Euler characteristic (and vice versa). This oscillation continues up to level 
E'^. It would be of interest to find a simple and intuitive geometric reason for this large 
oscillation. By symmetry above E'^ we have "essentially" covered the whole sphere minus 
an exponentially large number of disjoint simply connected sets. 
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Remark 4. The above Theorem also holds as stated in the pure p-spin case. Only the 
complexity function Q^iu) needs to be replaced by its analogue given in Theorem 2.8 of 
|ABC10j (see also Remark [5] below) . 



Proposition 4.1 and Theorem |4.2| are proven in Section [8} 

5. Complexity of critical points 

5.1. Main Identity. In this section, we introduce the main identity that relates the 
mean number of critical points of index k with the A;-th smallest eigenvalue of the Gaussian 



Orthogonal Ensemble. This identity, given in Proposition 5.1 , is the analogous of Theorem 
2.1 of |ABC10| and it is the first step of the proofs of Theorems 2.1 2.4 2.6 and 3.3 

We fix our notation for the Gaussian Orthogonal Ensemble (GOE). The GOE is a 
probability measure on the space of real symmetric matrices. Namely, it is the probability 
distribution of the N x N real symmetric random matrix , whose entries {Mij,i < j) 
are independent centered Gaussian random variables with variance 

MM^- = (5.1) 

We will denote by Egoe expectation under the GOE ensemble of size N x N. 
Let Xq < \^ < • • • < ordered eigenvalues of M^. 

Proposition 5.1. The following identity holds for all N , u, k ^ {0, . . . , — 1}, and for 

all Borel sets B cR, 



E[Crt^,,(S)] = C{N, u', u") 1^ E^aoE [ exp { y ((A^)^ " " ^ (A^ " } 



whereC{N,u',u") = 2^/^r^)^. 



(5.2) 



Proof. Proof of Proposition |5 . 1| is a rewrite of the proof of Theorem 2.1 of [ABClOj with 
one subtle difference: the law of the Hessian in the mixed case gains an independent 
Gaussian component on its diagonal. In this proof, we use H to denote Hj^ j^. 

The hypothesis on i' allows us to apply Rice's Formula, in the form of Lemma 3.1 of 
|ABC10j . It says that using dcr to denote the usual surface measure on S^^^{^/N), 

ECrtjvfc(5)= / E[\ det V'^H{a)\l{H{a) e NB,i{V^H{a))) = k}\VH{(T) = 0] 

(5.3) 
X (j)cr{0)da 

where (pa- is the density of the gradient vector of H. 

Now, since H is invariant under rotations, to compute the above expectation it is 
enough to study the joint distribution of {H, VH, V^H) at at the north pole n. We fix a 
orthogonal base for the tangent plane at the north pole, and we consider V H (n) ,V'^ H (n) 
with respect to that base. Denoting subscript by a derivative according to a orthonormal 
basis in TcjS^~^{^/N) we have that 

Lemma 5.2. For alll<i<j<N — 1, 

E[Hin)^] = N, E[H{n)H,{n)] = E[Hi{n)Hjk{n)] = 0, 

E[H{n)Hij{n)] = -u'6ij, E[Hi{n)Hj{n)] = u'Sij, 
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anc 



E[Hij{n)Hki{n)] = —[iy"{6ik6ji + dudjk) + {y" + 

Furthermore, under the conditional distribution ¥'[-\H{n) = x\ the random variables Hij{n) 
are Gaussian variables with 

X 



N 



and 



i.e., if M is distributed as a {N — 1) x {N — 1) GOE matrix 

E\V''H\H(n)] = (^^2i/")^/'M^-i + -^{aZ - -^u'H(n))I 
where Z is an independent standard Gaussian. 



Last Lemma implies that (5.3) can be rewritten as 
ECrt^,,,(5) 



E 



det ((^^^2^.")'/'^'^"' + ^(ViVaZ - u'H{n))l) 



^'^)l]=k\l{H{n)eNB] 



H(n) 



where ujn, the volume of the sphere S'^ ^{V^), and 4>n{n) are given by 

.1 2^^/2 



UJN = (VN) 



r(iv/2)' 



(5.4) 
(5.5) 



Since we can assume q / (the case a = 0, i.e. the pure p-spin was treated in |ABC10| ). 
we can rewrite the conditional expectation in (5.4) as 



Vn „iv-i. 



27r 



N 



JV-l 
2 



B 



e^E |det (M^-i - X{y))l\ l|i[M^-i - X{y)l] = fcjdy 

(5.6) 

where X{y) is a Gaussian random variable with mean m = (21/"^^^ 1^)1/2 ^^'^ variance 



= 2t//'(jv_i) • Hence, we can apply Lemma 3.3 of |ABC10j with G = M to get that (5.6) 
is equal to 



N ' 



'IT 



t2 



E, 



B 



N 

GOE 



, ^ ( ( \N\2 2 2;/" , 

expi w[(^k ) -y - -^\^k 



Putting (5.4), (5.5) and (5.7) together we end the proof of Proposition 5.1 



□ 



5.2. Proof of Theorems 2.1, 2.4 , 2.5 and 2.6 
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5.2.1. Proving Theorem 2.1 and Proposition^^ In this subsection, we will compute the 
logarithm asymptotics of the left-side of (5.2). 
Let F : - 



be given by 
1 / v" + v' 



2V2vV'zy' 



v" + v' 



f'?y ~^ II \ I 

ylA jjll _|_ yl 



,11 



\y 



v' + V 



12 



U" + v' 



,12 



Note that F{X, y) = —ay"^ + hyX — cA^ for some constants a,b,c> 0. Let 

^ \/z2 _ 2dz = - (^x\/x2 - 2 + log[2] - 2 log (^x + - 2 
For any /c G N fixed, let 



^llog'^ + FiXl[u],u)-ik + l)hi\Xl[u]\), if u<-E^ 
where ^^^^^Ig^ < Xl[u] < -^2 is given by 

""x2-2(A: + l)Ii(|x|), 



(5. 



(5.9) 



(5.10) 



-ux 



i.e., A^[n] is a solution on (— oo, — \/2] of 



z/V2i/" 



/2 



A^,M + (fc + l)./(A*M)2-2 = 0. 



(5.11) 



Our goal in this section is to prove that 9k,v is the /c-complexity function. When = 
the formula for Oq^i, simplifies as follows. 

Proposition 5.3. For all u G M, 



(n) 



V — v'-^+v 



u^(y'-2) 

4(!/'-l) 



^/2^v>{v'-\)' 



if - Eoo < u, 
if u< -Eoo- 

(5.12) 



Remark 5. It is possible to recover all complexity functions of the pure case by taking a 
to zero (i.e. recover the first results of |ABC10| ). In particular, if a = 0, = Eoo and 
we do not have the intermediate regions where the k- complexity functions are equal for 
different k and non-constant. 



We postpone the proof of Proposition 5.3 to the end of this subsection since we will 
need another characterization of Oi^. ^. 



Proof of Theorem 2.1. To prove Theorem 



rithm asymptotic limit of the left hand side of (5.2) 



2.1 



it suffices to show that 6k,uiu) is the loga- 



First, note that we can rewrite (5.2) as 



C;vEe-^^(^"'^-)l{Y;v G B}. 



(5.13) 



where is a Gaussian random variable of mean zero and variance N independent of A^ , 
E is the expectation with respect to GOE and Yiv and 



hm llogC7^ = ;^log4, A(A,y) = F(A,y) + ^ 



1 



A^ 



vy \2 

(5.14) 
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By the independence of Yn and and Theorem A.l of [ ABClOj . the sequence of 
random variables (A^, Yat) satisfies a large deviation principle of speed N and rate function 



4(A,x) 



^ + (A; + l)Ii(|A|), if X<-V2, 
oo, otherwise. 



Therefore, in view of (5.13) and (5.14), we can apply Laplace- Varadhan Lemma (see 
e.g. |DZ98j . Theorem 4.3.1 and Exercise 4.3.11) and get that 



lim llogECrt^,fc(i?) = J 



log — + max < A 

xeB,X<-V2 



-Xv'x-V2i7'Xf -2h{\,x) 



We will now analyse the above variational principle. We start by the case ol B = (— oo, u). 
We want to find 



max 

x<u,\<-'/2 



x^ + A^ 



1 



a" 



{u'x - y/27'Xf - 2{k + l)/i(|A|) \. 



(5.15) 



Case u > —E'^: If n > —E'^ then we maximize (5.15) in x first. The maximum is 



obtained at x = a;;^ := '^^/Y^^'' A < u. Plugging x\ back in (5.15), we get an increasing 



function in A, since /i(|A|) is itself decreasing. Thus the maxima is realized at 

x = x\, A = — \/2. 

—oo,u) with —E' < u. 



This together with (5.2.1) proves Theorem 2.1 in the case B 



Case u < —E'^: In the case u < —E'^, xx < u ii and only if A < ^w^. Therefore if x* 



maximizes (5.15) then 



xa A < 



\/2' 



u 



EL 



and X* =U'^ < A < -V2. 



(5.16) 



If we plug in the correspondent values of x in each region we note that in the first case 
our function is again increasing in A. Furthermore, since at A = ^w^, xx = u, we are led 



to the following variational principle valid in both cases of (5.16) 



max 

#i<A<-v^ 



-u^ + A2 - A,{i^'u - V2?^Xf - 2{k + 1)/(|A|) 



/2 



-(1 + ^)^2+ max <^ — uX ^ A2 -2(fc + l)/i(|A|) 



a" 



,'2 



(1 + -^)u^ + max ^'(A) = max r(A) 



\/2u 



f5±<A<-v^ 



(5.17) 



Note that ^'(A) is a parabola a}? + 6A,a < plus an increasing function. The critical 
point of the parabola is given by 

v'^f2A/'u 



X 



> -V2 



u > -En 



(5.18) 



Therefore if u > —Eoo, ^ is an increasing function in A, so its maximum is attained at 
A = -\/2. This proves the Theorem in the region —Eoo <u< —E'^- 
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li u < -Eoo, equation ( |5.18[ ) and the facts that ^''(-\/2) < and "^'{Xc) > imply 
that the maximum is taken in the interior of the interval [Ac,— v^] at ■^fei''^]- This ends 
the proof of the Theorem in the case B = {—oo,u). 

Now, it is easy to extend it to any open set B. Let u* be the point that realizes the 
supj^g^l 6k,uiu)- Prom the continuity and uniqueness of a local maxima of Of^^,^, it is clear 
that either u* = —E'^ or u* is in the boundary of B. Assume without loss of generality 
that there exists an increasing sequence Un in B approaching u* . Since B is open, there 
exist en > such that 

E(CrtAr,fc(-00, Un) - CrtAr_fc(-00, Un - Cn)) = E CTtN,k{Un " Un) < E CrtN,k{B) 

< ¥,Crt]\f^k{—oo,u*). 
But since 6k,u is continuous and increasing for u < —E'^ last equation implies 



Ok,u{Un) < lim —logECTtNk{B)<dk,u{u* 
Af— s>oo iV 



for all n, which proves Theorem 2.1 for any B open. 



□ 



It remains to prove Proposition 5.3 We first need the following miraculous Lemma. 



Lemma 5.4. For all u < —Er, 



d 



Proof. The proof relies on how we derived 6o^i,{u). When u < —E^o, Oq^u) is the maxi- 
mum over A of a functional P (that depends on v") given in (5.17). Its maximizer \*{u) 
is the smallest root of a second degree polynomial that can be derived from (5.11). This 
second degree equation is given by ^4 + BX + CX^ 

2u'^v'^v" 

A = 2 + - 



where 



B 



C 



2^/2uv' yJV' {{-I + v')v' + u") 
2 ((-1 + v')'^u'^ + 



(5.19) 



Now chain rule and the fact that A*(ii) is a maximizer imply that ^ 



only if gf,(r(A*(n)) 



if and only if 



P (A*(n)) 



g^b'o,;/(n) = if and 
0. The last condition can be 



written down as a second degree equation of the form 




u 



v") 



u 



+ 



U' + U'^ + U") A2' 



+ 



ly'^ + 1^" {u' - i/'2 + ly") 



2y/2uv'^X 



0. 



+ 



^/2uv'X 



(5.20) 



Comparing the coefficients of (5.20) with (5.19) one sees that their ratios are constant 
equal to This immediately implies that they share the same roots. So A*(ti) indeed 



satisfies ^g^T/P^ {X*{u)) = and the lemma is proven. 



□ 
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Proof of Proposition 5.3 From Lemma 5.4 we know that for u < —Eoo, Ok,u does not 
depend on v" . By choosing v" = v'"^ — u' -\-e and taking e to zero we get the desired result. 
Indeed, when e goes to zero 



A*(n) 



uu 



, F(X(u),u) 



-u^(v' -2) 
- 1) 



□ 



5.2.2. Proof of Theorem 2.4 We want to prove that there are no critical values of index 
k of the Hamiltonian above —N{E^ — e). The function 9k^u is strictly decreasing on 
{—E^,oo). Using Theorem 2.1, we have 



IE [ Crtiv,fc(-^- + e, oo)] < exp {Nek,,{-E^ + e) + o{N)}. 

The constant —E^ is defined by 9k,v{—E^) = for all k. Therefore, 9k,p{—Ek + e) 
c{k, ly, e) < 0. An application of Markov's inequality as 



B^,,(e)) <E[Crtjv,fc(-i?oo + e,oo)] < e'^'^^^''^'^) 



proves Theorem 2.4 for the event BN,ki^)- The proof for the event ^iv,fc(e) is analogous. 



5.2.3. Proof of Theorem\2.^ The proof of Theorem 2.5 follows the same steps as the proof 



of Theorem 2.1 First by Lemma 3.5 of [ABClOj . for any e > 0, there exists a constant 
c = c(7, e) > such that 



lA 



N 



>e)< 



Therefore if we use Proposition 5.1 (5.14) and the above statement we have that for 



any e > 0,6 > there exists constants c = c(e), d = d{e) such that for N large enough 

JB 



JB 

< CNe^ ™PAe(s^-.,s^+.),aei3{^(A.s/)}(]_ + 6) + e'^^e"''^^ 



On the other hand we have the lower bound 

Jb 

inf 



Crt^,fc(i?) >Cn I e"^ {n^^^y)) i{Af e {s^ - e, + e)} 
Jb 

,_£,s^+E){-F(A,y)}^^ 

JB 



Taking log on both bounds and taking e to zero afterwards, we see that 



iV 



logECrt^,fc(S) = sup{F(s^,y)}. 



yeB 
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5.2.4. Proof of Theorem 2.6. We now prove the asymptotic limit of the mean number of 
critical points at some level of energy. 

Since the total number of critical points is greater than the number of critical points of 



index k{N) with k{N) satisfying (2.21) for 7 G [0, 1] we clearly have the lower bound 

(5.21) 



1 



sup sup 6'-y ^(ti) < lim — logE CrtAr(-B). 

7e[o,i]«eB ' N^oo N 

For u < —E'^, taking 7 = (i.e. considering the complexity of local minima) we get the 



right hand side of (2.24). For u S {—E'^,E'^) the supremum on 7 of 9^^u{u) is attained 



at 7 G (0, 1) such that = ^w^, plugging it back on the left hand side of (5.21), we get 



the right hand side of (2.24). Last, for u > Eoo, one just need to take the complexity of 



local maxima. This is enough to prove a lower bound. 

To show a matching upper bound, we proceed as follows. A sum over k in Proposition 



5.1 , gives us that 



E[Crt;v(i?)] = 2iVy^(^)f ^E^OE / exp |7VF(z, y)|dyL 



N{dz). 



and -Ljv is the empirical spectral measure of the GOE matrix. The constant in front the 
integral gives a constant term C^, after the log limit. Furthermore, 



/ E^oE f exp\NF{z,y)]dyLN{dz)<N [ supexp |iVF(z, y) Idy 
Jb J I ) Jb zm L J 



AT u" 



(5.22) 



IB 



Soii Bn{-E'^,E'^) ^ t his matches the right hand side of ( |2.24D . If B O {-00, -E'^ 
then we can estimate ( |5.22 ) with 



N 



B 



^GOE 



/exp{ 



NFiXo,y)\. 



Applying log, dividing by and taking limits we get Theorem 2.6 from Theorem 2.1 



6. Partition Function 



In this section we prove Theorem 3.1 and we derive a formula for the 1-RSB solution at 
zero temperature that will be useful in the next section. We start by the proof of Theorem 

EH 

Proof of (a): By Holder's inequality the function jr'Elog Z]\[^,j{(3) is convex in (3, therefore 



its limit -Foo(/3) is also convex. From (3.6) 



< liminf 4^00 (/3) < hmsup < ^o- 

/3-s>oo p ^_^oo P 

So F{l3) is convex, positive and grows at most linearly. This easily implies that 

lim ^Foo(/3) = sup^Foo(/3) G [0,oo). 



(6.1) 

(6.2) 
□ 
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To prove item (b), we will need to introduce some notation and the proposition below. 
Let <T* be a point on the sphere such that Hj^^i3{a*) = NGSjy and let d denote the 
geodesic distance on the sphere. For p,a, K > 0, let 



Bn,p = <; ct g Sn-i{VN) : d{(T,cT*) < p 
and A^^a,K{N)^ be the event 



sup |i/jv,/3(c^) - NGSn\ < KNe"" 



(6.3) 



Lemma 6.1. For any < a < 1 there exist constants K, Ki > so that for all e > and 

all N sufficiently large 

'A,a,K{Nr) < 2e-^^^. (6.4) 



Note that this bound is independent of e. 
Proof. Clearly, 

A,,KiN) 5 Ao^^KiN) = |||/7jV,/3||a < KN^-': 

where 



\HN,l3\\a = sup 



(6.5) 



by 



d{a,a'Y 

Now consider the centered Gaussian process Xa held on Sn-i{VN) x SN-i{-\fN) given 



0, otherwise. 



Since the Gaussian field H^^f^ is almost surely, then 



But now a simple computation yields for a ^ cr\ 



2N 



d{ai,a^ 



I \2a 



2N 



NO) 



2a 



1 — z^(cos 9) 



(6.6) 



(6.7) 



(6.8) 



where 6 is the angle between <t,<t' in M^. 

Therefore by the boundedness of z/'(x) in [—1, 1] there exists a constant C independent 
of such that (if a < 1/2 or a < 1 - using the boundedness of z^'(x) and i'"{x)) 

sup EXl{a,a') < CN^''^. (6.9) 



Now, by Borell's inequality, (see page 50 and 51 of |AT07j . where we take u = KN 
o'T ^ CN^~°') for all 5, if N, K is large enough 



M supX„(<T,cr') > i^iV^-t ) < e 

V <t,<t' 



.10) 



Taking i^i = K"^ /AC in the last equation, using (6.7) and symmetry of Xa the lemma is 
proven. □ 
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Now we can start the 

Proof of (h): We will show that for any 5 > there exists e{5) so that if N is large enough 

(6.11) 



The proof of (b) will then follow from (6.11 ) and Borel-Cantelli's Lemma since for all (5 > 
by Lemma|6.1| 



Y,¥(\GSN + f^\>6 \ < 

N=l ^ ^ 



oo. 



(6.12) 



We will prove (6.11) by showing that for any > if is large enough A^^ct,K{N) C 



{\GSN + foo\<S}. 

On A,^a,KiN), 

ZnA^)= f e-/^^-./'(-)A^(d^)>e-'^^«^--^^^^"A^(i?^^J. (6.13) 



Recall that A7v(d<T) is the surface measure of SNivN) normalized to be a probability 
measure. We trivially have the bound 



1 

N 



log ZnAI^) < -PGSn. 



(6.14) 



Combining (6.13) and (6.14) we then have on A^^a^xi^), 



-^logZ^AP) - Ke'^ + -i_logA^(i?^^^J < GSn < - log Z^,,(/3). (6.15) 

Note that using spherical coordinates and the inequality ^ < sin0 for 6 < |, we have 
for e < tt/2, 



sin^-2(0) ^<Pji^J sm'^^\^) dcj) 



> 



N-l 



1 



(6.16) 



•vr' 7r(iV-l) 
So on A^^a,K{N), for some constant C > 
1 



1 



log ZnAI^) - i^e" + Ce < G5iv < log Z^,.(/3). 



Therefore by (3.2), for any > one can take N large enough so that, 

i^oo(/3) 



/5 



f3 P P 



.17) 



(6.18) 



By taking (3 large enough, part (a) of this Theorem and by choosing e sufficiently small, 
(6.11 ) is proven. □ 



7. Proofs of Theorems 3.4, 3.6 



AND 



3.8 



In this section, we prove Theorems 3.4 3.6 and 3.8 
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7.1. Calculating /i. We now prove equation (3 
Lemma 7.1. 



1) There exists e, M > such that the following limit holds 



fi 



lim iFi(/3) 



i„f 1 

a,b(^[e,M] 2 



r / 1 /, « + ft 

6 + a + 7 (log 

a 



(7.1) 



(2) fi depends (continuously) only on the first derivative v' . 

Remark 6. It is remarkable that while the fc-complexity function depends on the first two 
derivatives at 1 of the covariance function v and /i depends only on the first derivative v' 
and Eq{v) = fi for any pure- like mixture. 

Proof. First, taking fi = m6r + (1 — m)5q, we can write the IRSB Free Energy Fi{P) via 
the Crisanti-Sommers representation (see (1.4) of |TP07j ) as 



T 1 

1 TT \ ( log(m(g 

m,r,qe[0,i]3 [ m{q — r) + 1 — q m 



inf 



r) + 1 



log(l - q)) + log(l -q) + p^m{u{q) - u{r)) + p\v{l) - u{q))\. 



It is easy to show that the infimum above is attained at r = 0. Therefore, 



inf 

m,,ge[0,l]2 



— log(- 
m 1 



1 



q + mq' 



(7.2) 



The conditions to be a critical points are 

p-'v'iq) = 

. ' 1 - q + qm 

log 



(1 - q){l - q + qm) 

2a2 . 



u{q)m^fi'^ + 



(7.3) 



1 — g / 1 — g(l — m) 

Let {q*,m*) = {q* ,m*){j3) be a solution of \7.J>\ . First, from the first equation we deduce 
that either q* goes to 0, 1 or 1 — 5*(1 — m*) goes to zero as /3 goes to infinity. Analogously 
to the pure case, the case q* approaching zero is excluded via Lemma 3 of |TP07j . Since 
1 — q*{l — m*) — )• implies g* — )• 1, it is then a fact that q* ^ 1. Looking at the second 
equation we see that m* has to go to zero. Indeed, if it is not the case 1 — ~ 
implying log ~ which is not possible. 

If we put A = 1 — q* , B = m* , then A, B go to zero as /3 goes to infinity and conditions 
(7.3) become 



1 



A^ + AB' 



log 



1 + 



B 
~A 



B^P^ + ( 1 + 



B 
A 



-1 



(7.4) 



We will argue that B A, i.e. lim^_^oo ^ = ^ G (0, oo). Suppose that A >> B. The first 
condition in (7.4) implies B « But now the RHS in the second condition goes to 

1 while the LHS goes to 0. Next, suppose that B >> A. First condition in (7.4) implies 
— while the second implies log(l + ^) ~ ^ which is a contradiction since ^ 



AB 



oo. 



Therefore, not only we have B A but we can also write 

A = a{(5)a* B = h{P)b*r^ (7.5) 

where a(/3), 6(/3) are functions that converge to 1 as /3 goes to infinity and a* ,h* G (0, oo). 
Furthermore, a* ,h* satisfy the following relations: 



1 



a* (a* + h*] 



log 



1 + 



h* 



ib*f + ( 1 + 



b* 



(7.6) 
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To get the statement just note that replacing m = bP and g = 1 — a/3 in (7.2), we 



get 



where 



P{/3,a,b) 



+ 



(/3 - 6-1) log (1 - (1 - a/3-1)) ^ [1 _ (1 _ a/3-i)(i _ 5/3-^)] }. 



Clearly, for any a, 6 > the function P(/3, a, d) converges pointwise as /3 goes to infinity 
to h(b + y'a + ^(lo| 



a+b ' 
a ' 



|. Since we know that the location of the minima of ^-Fi(/3) 

converges to (a*, 6*) € (0,oo)^, this is enough to guarantee the convergence stated in 
Lemmma 17.11 



By solving for the critical points of (7.1), i.e. using equations (7.6) we can get an 
expression for /i in terms of v' . Namely, 



/l 



v'y 



1 1 

u'y y v'y"^ — 1 



log(i/'y^) ] = y + 



1 



where y = y{v') is given by the unique solution of 



u'y"^ - 1 \ 2 
v'y 



y + 



v'y' 



v'y 



ylog(z^V) 



In other words, y = ^= where a is the unique solution of 

alog[a] - a + 1 - ^ ^ = 0, 
This ends the proof of Lemma |7.1[ 
7.2. Proof of Proposition [3T3[ 



y>v 



a> 1. 



yv' 



r-l/2 



□ 



Proof. If V is pure-like then 9k,v{—Eoo) > 0. Since 6k,u{u) converges to negative infinity 
as u goes to negative infinity, Ek{v) are well-defined. Furthermore, as k goes to infinity, 
A^(ii) converges to — \/2 for any u < —E^o, implying that 9k^y{u) converges to F(-\/2, u) 
pointwise. Therefore, taking u in a small neighborhood of E^ and using the fact that 6k,v 
are increasing in that neighborhood we see that the zero of 9^^^ has to converge to the 
zero of F(-V^, n). Namely Ek (u) converges to E^. □ 

7.3. Proofs of Theorems 3A[ |3.6| and 3.8 We start by the case when v is critical, i.e. 
by the case when G{v' , v") = 0. 

Proposition 7.2. A mixture v is critical if and only if 

-v' + v'^ 



fl — Eoo — Eq^u 



(7.7) 



Proof. If V is critical then y = ^^-^ is the unique solution of (7.1) with y > -1=^. Indeed, 



" {-V' + V") {-V' + + u" 



+ 



-1 + ^ 



/3 



0. 
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Plugging back the value of y in (7.1 ) we get fi. On the other hand, if / 



then 



one solve equation (7.1) in y to see that the only positive solution is y = By the 

definition of y in (7.1) this immediately implies that is critical. And trivially, i' critical 
is precisely the case where E^o = Eq^^. □ 

Now we analyse the case where is critical or a full mixture, i.e. the case where 
G{i'' , v") < 0. In this case, the zero of the complexity function can be explicitly computed 
and is given by: 



where E^ was defined in (2.14). Note that Eq^ is a function of u' and i'" . 



Proposition 7.3. IfG{u',v") < then 

-JiEq^^ = if and only if G(z^', v") = 0. 



Proof. Let 



Calculating the derivative -J^Eq^u one gets 



A{u', u") = Mu" - + v') {u' + v") lo, 



2{v" - v') 



u V {v + I ioe 



X ( 2v"(v' + v"fA{y\ v") 



-1 



(7.8) 

Sufficiency comes from a simplification of the above formula. To get necessity we solve 



a second degree equation on the variable M = log 
equation has a unique zero given by 



to see that this second degree 



This is precisely u") = 0. 

With the above propositions we now prove Theorem |3.8| 



□ 



Proof of Theorem 3.4- If is critical Theorem |3.8| is Proposition |7.2[ Now suppose that 



u is pure-like. By Lemma 7.1 and (5.10), both /i(i^) and Eq{v) are independent of u" . 
Consider then another mixture /x such that /i' = v' and ^ satisfies /i") = 0. Since G 
is continuous on its domain, we have 



f^{u) = f^{^j.) = E^{li) = E^{v). 



□ 



Proof of Theorem 3.6. The proof is just a simple but major simplification. Note that by 



solving (3.12) in a we can rewrite (3.12) as 

2Log 



1 



-h{-2 + v)- + 



i ( 2 + b^i^' - bV^V^+Wi^' 
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This implies that a critical point b*(u) of M{u, b) = ub — bfi{b) is given by 



6i(n) 



-u{i'' — 



2{v' - 1) 

Here, we choose b\(u) since b*_{u) is negative and a local maxima. Furthermore, since 



(^M)(n,c, 



n + 2 



1 



<0, 



b\{u) is a global minima in [c(z^),oo). So, replacing in (3.12) and setting z 

—u(v' — 2) + ^fu'y/ v?v' — Au' + 4 = we get the function 



zu 



+ 



2(i/' - 1) ' 8{iy' - 1) I 2(z^' - 1) 
1 / 



v'z^ 



4(1/' - 1)2 



i^'z^ 



+ 



V Z 



4(z/'-l)2 2(i/'-l) 



'4 + 



A{v' - 1)2 



(7.9) 



Comparing term by term the logarithm, polynomial and fractional terms of (5.10) and 



(7.9) one gets the first part of the Theorem. The proof of the second part follows from 
Fenchel's duality theorem (see |Roc70l Chapter 1]). □ 



Proof of Theorem 3.8, First, note that for a fixed u' , there exists u"'^ such that the condi- 



tion that v is a full mixture can be written as i^" > v"^. From Proposition 7.3 we know 
that 3^77 -Eo.i/ does not change sign as we change u" . Taking i/" to infinity in (7.8) we see 



that gp/E'o,^ > 0, meaning that Eq^i^ is increasing in . Since at u'"^, £'o,i/ = /i; we get 
that for any u such that G{u' , u") < 0, £'o,v > fi proving the Theorem. □ 

8. Proof of Propositon 14.11 and Theorem 14.21 



In this section we prove Propositon 4.1 and Theorem 4.2 



Proof of Proposition 4-1 We start from the following identity: 



E x(.Au) = Y.{-l)'CTtk(,A{u)) 

X^(-l)''' / lEM det V^i?iv(cT)|l{j(v2Hjv(fT))=fc}l{HAr(a)<7V«} 



X <^v//^(0)dc7 = {2u'7r)-'^\S''-\VN)\ 



VHn{(j) = 



^ / E (-l)VetV2i?^(^7)|l|,(v2^,^(.))=,} 

i._n J-oo \ 

J — OO \ 



1 2 

Hn[(^) = X )e~2iv^ dx 



, , iv-i 2tt 2 iv-i VA'" 



N 2 

HN{(y) = Nx e-T^ dx. 
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Lemma 8.1. If Mn is a N x N GOE with variance EMfj 



2N 



then for any x E 



Edet(M^ - xl) = 2-^iV-t(-l)^U^(ViVx) 



iV 



where h]\f(x) is given in (4.2). 

Proof. The proof, a straight-forward Unear algebra exercise, can be found as Corollary 
11.6.3 in [^07]. □ 



Now by Lemma 5.2 



JV-l _ J V-2 A'" 2 



27r 



OO fU 



N -I 



N 



El det 

-OO J — OO \ 

The double integral becomes: 

coo fU / 



2u"f'^M^-^ + {ay - u'x)I 



■1) 



_N 2 _ N 2 

e e dxdy. 



^ N ' 



— OO J — OO 



Ej^det 



M 



N-l 



+ ( 



>jV- 1 

' N 



2u")-2{ay-u'x)I 



N 2 N 2 

e~T^ e~^y dxdy, 



which by Lemma 8.1 can be rewritten as 



~ , ,VN{u'x-ay). _N^2 

hN^ii ^F^^ -)e 2^e 2 



A[^2 



OO ./ — OO 



'2u" 



dxdy. 



Combining (8.1) and (8.2) we get Proposition 4.1 



□ 



We will need the following Lemma to prove Theorem 4.2 



Lemma 8.2. Let a, h he constants such that a > 1/2 and 6 > 0. Set 



In{M) 



M 



(/.jv-i(\/iVa;)e'^^'+^^ dx. 



As N goes to infinity: 

(1) IfV2<M then In{M) = 0(e-JV(aA/2+6Af+/i(Af)))^ 

(2) // -\/2 < M < \/2 and if we set M = \/2cos w with e < uj < ir - e then In{M) is 
equal to 

2l/4 g-7V(aM2+6M) 



sm 



7ri/2iv! 2|m'(2i(M))|(sinw)i/2 
(3) IfM<-V2 then In{M) 



■ N 1,, , 37r , , 

- - -) ( sin 2l^ - 2^) + — + a{M) 



(l + 0(iV-i)). 



_|_g NX{a,b.M) y;/jgj^g X{a,b,M) is the minimum of 

ax^ + bx + Ii{—x) in [M, —y/2] and c is a positive constant that depends on a,b 
and M. 

A few comments before the proof of the above lemma. First, under the assumption 
that a > 1/2 and 6 > the major contribution to the integral in part (2) comes from 
a small neighborhood of M, instead of the minimum of ax^ + bx. This is due to rapid 
oscillations of (j)N-i inside the "bulk" {—^/2,\/2). Second, in part (3), the condition that 
the minimizer of ax^ + bx + x) lies inside [M, — \/2] is similar to the condition on 
(5.11). It will lead to the asymptotic Euler characteristic in the region u < —E'^. t The 
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main tool to prove Lemma 8.2 is the following well-known formula for the asymptotics of 
the Hermite functions, first proved by Plancherel-Rotach ^PR29j . Let 



h{x) 



V2 



V2 



1/4 



+ 



+ \/2 



V2 



1/4 



Lemma 8.3 (Placherel - Rotach). There exists 5o > such that for any < S < 6q the 
following asymptotics hold uniformly in each region: 

(1) Ifx < -V2-5, 



Nh(-x) 



h{x){l + 0{N-^)). 



(2) If-V2-6 <x< -V2 + 6, 



1 



,Af-l 



X + V2 



X-V2 
X + V2 



1/4 



X 



V2 



Ai (1 + O(iV-i)) 



'^'|^/i(-x)r^ Ai' (1 + O(iV-i)) 



where Ai(x) is the Airy function of first kind, Ai(x) = ^ J^^o (y ~'~ and 

(3) // —^/2 + 5 < X < \/2 — 6 and if we set x = Vlcosu with e < lo < tt — e then 

ol/4 1 / N ^ 37r\ 

<^n-i{VNx) = — ^- ^ sin (- + )(sin2^ _ 2^) + (1 + 0{N-')). 

TT^/^Ni (sina;)2 V 2 4 A J 

(4) IfIfx>V2 + 6, 



-NIi{x) 



h{x){l + OiN-^)). 



Proof of Lemma\8.S\ Part (1): We can use the uniform asymptotics given by the ex- 



ponential region (4) in Lemma 8.3 . Precisely, by hypothesis, the function K{x) :- 
ax"^ + hx + Ii{x) is increasing in [M, oo) and by Laplace's method: 



In{M) 



oo ~N{ax^+bx+Ii{x)) 



M \/47r\/2iV 

^-NK(M) 



h{x){l + 0{N-^))dx 



N\K' {M)\^/aW2N 



h{M){l + 0{N 



Part (2): Choose 5 < Sq such that -^/2 < M < y/2 - 5. We split the integral In{M) 
into three parts 



In{M) 



+ 



M 



+ 



V2-5 



■.= h{M)+h + h. 



.3) 



We will show that the main contribution in this case comes from /i(M). As in part (1), 
it is easy to see that 

J3 = 0(e-^^(v^)). (8.4) 



COMPLEXITY OF RANDOM SMOOTH FUNCTIONS OF MANY VARIABLES 



30 



Next since \x |i/^|Ai(x)| and |x|-i/4|Ai'(x)| are bounded functions on M a change of vari- 
ables z = Ii{—x) when using Part (2) of Lemma 8.3 immediately implies that for any 
e > 0: 

I, = 0(e-^(a(V2-5)2+6(v^-5)+.). (8.5) 

Now we estimate Ii(M). Using the uniform asymptotics of i;^>Ar_i we need to evaluate 



2i/4 



i/SATj Jm 



V2-5 



-N{ax'^+bx) _ 



1 



(sino;) 



Y sin ( ( T)(sin 2u: — 2(x>) H — — ] dx. 



iV _ 1, 
2 4' 



4 7 



.6) 



Performing the change of variables x = -v/2costj,0 < uj < tt the integral above becomes 
(for some different S > 0) 



-N{2a cos^ u}+V^b cos oj) 



(sincj) 2 sin ( ( 



.N 



1 Svr^ 
-)(sin2a; - 2a;) + — jdw (8.7) 



for l{M) = arccos(2^/^M). We now rewrite cos^w = i+^o'^ and use the substitution 
2a; = 2 to obtain the integral 



1 

71 



2n-2S 



I2l(M) 

Last, we write 

,N 



-Ar(a+acos2+4=cosf). . Z l . ( N 1 SvtX 

v2 2/(-gjj^ _j2 sm ( — -)(smz — z) + — Idz. 



sm((y-^)(sin.-.) + ^)=l 



, (8.9) 



where /i(2) = —\{smz — z) + ^ . 

Therefore, we just need to evaluate the asymptotics of 



2n~25 



2l{M) 



-Nm(z) 



j{z)dz, 



2t(M) 



-Nn{z) 



k{z)dz 



where m and n are entire functions given by 

m{z) = a + acosz-ir 



b z i 

cos - — -(sm z — z) 



^/2"""2 2' 

h z i . . 

nix) = a + acoszH ^ cos — I — (smz — z 

^2 2 2' 



(8.10) 

(8.11) 
(8.12) 



and j{z) = sin(|)5e^-^i(^), k{z) = sm{^)h-^^^^^l 

We will change our contour of integration and apply Laplace's Integral in the appropriate 
integrals. Notice that the steepest descent paths are given by the equations 

■ ( ■ ^ coshy\ b X y X 
^•(771(2;)) = smxIasmhyH IH — ^ sm — smh = constant (8.13) 

V 2 / y2 2 2 2 

f coshy\ b X y X 
ks[n[z)) = smxl asmhy JH — ^ sm — smh — | — = constant. (8.14) 

V 2 / y2 2 2 2 

The phase diagram for the steepest paths of m is described as follows. First all lines 
X = 2k'K, S N are steepest paths. Second, for every t G (0, 27r) the steepest path that 
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passes through t goes from — ioo to vr + too if 6 > and from vr — ioo to vr + ioo if 6 = 0. 
The real part of m(z) is given by 

(1 \ b X y 

a cosh y H — sinh y \ + a -\ — — cos — cosh y (8.15) 
2 J ^2 2 2 

/ 1 \ b X y 

K(n(z)) = cos a; I a cosh y sinhy)+aH — ^ cos — cosh yH — . (8.16) 



2 y y/2 2 " 2' 

If we integrate m{z) between two points a,/3 G (0, 27r), we can deform our contour to 
be equal to the two steepest paths that connect a and (3 to z = — ioo. Precisely, we 
deform our countour into three pieces: we first follow the steepest descent path from a to 
a point with imaginary part yg < 0, |yo| large. From there we go along the straight line 
y = yo until we reach the steepest path that passes through /3, jy^, and then we integrate 



on this steepest path back to /3. From (8.15) we see that if we choose |yo| large enough, 



every point in the straight segment y = yo that we cross has real part x sufficiently close 
to so cosx > 0. This together with a > 1/2 implies that 5R(m(z)) diverges to infinity as 
y goes to negative infinity. The trivial bound 

e-^'^^'^j{z)dz\ < j e-^^(™("))dz sup \j{z)\ (8.17) 

combined with the bounded length of show that the contribution of this part can be 
made as small as we want by choosing yo large enough. 

In the two remaining paths the imaginary part of m is constant and therefore we can 
apply Laplace's method to get the asymptotic behavior. Since we assumed that M < \f2 
the contribution at 27r — 25 is negligible compared to the one at 2t(M). Indeed, by formula 
(7.2.11) of [BH9T], 

,-2-K-2& -Nmi2i.{M))+ii7T-aiM))-(2,(M)) 

/ e-^-(^)j(.)dz = + ))' (^-18) 

J2,(M) N\m'{2L{M))\ 

where a{M) is the angle of the steepest descent path of m at z = 2t(M): 

/ , / 1 — a cos z \ , , 

a{M) = arctan ( — j—^ ) . (8.19) 



2a sin z + 



V2 

The above argument adapted to the function n implies 
J2i{M) N\n'(2i{M))\ 



Noting that for any x G (0, 27r) |n'(x)| = |m'(x)|, we can combine (8.9), (8.18) and 



(8.20) to recover that Ii{M) is asymptotic equivalent to 

\l + 0{N-')). 



2I/4 ^-N{aAP+bM) 
Sill 

;ri/2ivl 2|m'(2i(M))|(sinw)i/2 



■ N 1., .Sir , , 

- - -) ( sin 2a; - 2a;) + — + a(M) 



(8.21) 

This ends the proof of Part (2) of Lemma. The proof of part (3) follows from the proof of 
part (2) and Lapace's method as in part (1) applied to the integral 

'—V2-S 

e"^'+^^+^i(-^)/i(x)dx = 0(e-^^(*^'"''')). 



'M 

We leave the details to the reader. □ 
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We now turn to the proof of Theorem 4.2 



Proof of Theorem We can rewrite (4.4) as 



oo J —oo 



where 



c(iV, ly) = 2iy"{[N 



JV-l 

2 — 

For the case a 7^ 0, we can change variables z = u'x — ay, w = ax + i^'y to get 

1 \ r , \/ 1 

u z + aw 



.22) 



Q,2 _|_ j^/2 



y = \v w — az 



and the above double integral becomes (using a = u" + u' — u ): 



1 



iVi/"(z^+m^) 



i/'z+aw<(u"+iy')^lL=ff 



So we have to evaluate the asymptotic behavior of the following integral: 



J= (PN-i[VNz)e 

'—00 .^—00 

r-M , TOO 



e dit; dz. 



We write the outside integral dz as /j^^^ + with M = ^^y^J," • The inside integral 
is just a Gaussian integral and therefore after a straight-forward computation the problem 
amounts to compute the asymptotics of the three following one- dimensional integrals: 





f'OO 


Jl = 






Jm 




POD 


J2 = 






Jm 



(l)N^i{VNz)e dz 
as J= (Ji + J2)(1 + 0(A^~^/^)) if iVis even and J = ( Ji - J2)(l + 0(iV"^/^)) if iV is odd. 



Take n < 0. We use Lemma 

v^2i^i/'m 



and 5 = - )^r7"2 > 



8.2 



in both cases. Note that by (|2.9p, o 



2(,y" + ;/'-jy'2) 



Now the condition M < — \/2 (M > — \/2) is exactly the condition u < —E'^ (n > 
—E'^). Applying the appropriate cases of Lemma 8.2 we see that the integral J2 is 
negligible compared to Ji. A comparison with ( |2.23 ) and (5.17) gives the proof of part 
(1) and part (2) of the Theorem with a and b as above, 



a(w) = arctan 



1 — a cos oj 



2asinc<j + 



b sin(a;/2) 
V2 



and f{io) = ( |m'(2a;)| sin^/^cj 



(8.23) 



where m is given in (8.11). Part (3) follows from symmetry of the Hamiltonian and 
(|431). □ 
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